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Abstract: We report on a complete study of the single top-quark production by direct 
^ . supersymmetric flavor-changing neutral-current (FCNC) processes at the LHC. The total 

cross section, o"(pp(gg) — > tc-|-tc), is computed at the 1-loop order within the unconstrained 
ly-^ , Minimal Supersymmetric Standard Model (MSSM). The present study extends the results 

I of the supersymmetric strong effects (SUSY-QCD), which were advanced by some of us 

' in a previous work, and includes the computation of the full supersymmetric electroweak 

' corrections (SUSY-EW). Our analysis of (T(pp(gg) — > tc -|- tc) in the MSSM has been 

■ performed in correspondence with the stringent low-energy constraints from b sj. In 

the most favorable scenarios, the SUSY-QCD contribution can give rise to production 
^ , rates of around 10^ events per 100 fb"^ of integrated luminosity. Furthermore, we show 

that there exist regions of the MSSM parameter space where the SUSY-EW correction 
becomes sizeable. This could be important, especially if the SUSY-QCD effects would be 
suppressed. In the SUSY-EW favored regions, one obtains lower, but still appreciable, event 
production rates that can reach the 10^ level for the same range of integrated luminosity. 
We study also the possible reduction in the maximum event rate obtained from the full 
MSSM contribution if we additionally include the constraints from ~ I^s- However, 
we treat these restrictions at a different level from the b — > S7 ones, due to the higher 
uncertainties inherent in the calculation of the matrix element associated to that mixing. 
In view of the fact that the FCNC production of heavy quark pairs of different flavors, 
such as tc or tc, is extremely suppressed in the SM, the detection of a significant number 
of these events could lead to evidence of new physics - of likely supersymmetric origin. 
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1. Introduction 

The forthcoming generation of high energy colUders, headed by the Large Hadron CoUider 
(LHC) at CERN, and followed by the future linear collider, depicts an exciting scenario 
for probing the existence of physics beyond the Standard Model (SM) of strong and elec- 
troweak interactions [1]. Among the possible discoveries envisioned for the physics at the 
LHC (some of them of a rather exotic nature, such as extra dimensions [2] and black-hole 
production [3]), we have the possible confirmation of the fundamental Higgs mechanism 
of Electroweak Symmetry Breaking. This would be accomplished in practice through the 
physical production of one or more Higgs boson particles. Undoubtedly, the next-to-most 
important discovery expected at the LHC is the finding of supersymmetric particles. 

Actually, the discovery of Supcrsymmetry (SUSY) (sec [4] for a comprehensive review) 
is intimately connected to the structure of the Higgs mechanism. In fact, unearthing 
supersymmetric particles would be strong evidence that Higgs bosons (in plural) should be 
around the corner. The opposite, however, is not necessarily true, but if a light Higgs boson 
of, say, 130 GeV would be found at the LHC, the hopes for SUSY physics would stay high 
and we would immediately felt encouraged to search for more Higgs bosons and potential 
supersymmetric particles. It is well-known that a light Higgs boson (m/j < 140 GeV) is a 
trademark prediction, if not of SUSY in general, at least of the Minimal Supersymmetric 
Standard Model (MSSM) in particular, which is after all the canonical scenario for low- 
energy SUSY phenomenology [5-7] . 

If SUSY is realized at the TeV scale (usually taJien as the characteristic energy scale 
to explain the naturalness problem of the SM [6]), one expects that a few (or even a 
bunch of) supersymmetric particles of the MSSM spectrum should be well reachable at 
the LHC. However, the tagging of heavy new particles is not an easy task because of the 
many decay modes available, most of them carrying invisible neutral species (some of them 
also of genuine SUSY origin, like sneutrinos and neutralinos) and, therefore, leading to 
missing energy events - usually hard to interpret. For this reason, one expects to get a 
complementary clue to the underlying SUSY dynamics from the short-distance quantum 
corrections on more conventional processes. If these supersymmetric quantum effects can 
be measured, they can be a solid handle to the properties of the new physics. The idea 
has been known for a long time and has been applied to the famiUar physics of the W and 
Z gauge bosons, see e.g. [8-11]. Here we wish to apply this method to the realm of rare 
processes, namely processes with conventional initial and final states which, although not 
strictly forbidden, turn out to be highly suppressed within the SM context. Among them, 
we have the fruitful Flavor-Changing Neutral-Current (FCNC) processes. 

The study of the flavor-changing interactions, in particular the FCNC processes, has 
been a very active field of research for about forty years, namely as of the glorious times 
when Glashow, Iliopoulos and Maiani (GIM) successfully proposed the existence of a fourth 
species of quark, the c-quark, to suppress to an acceptable level the strangeness-changing 
neutral-current effects in rare processes (e.g. fj,'^fi~) that otherwise would proceed at 

the tree-level, and similarly to further suppress the one-loop contributions in e.g. the K'' — 
system. Indeed, it was the experimental evidence that the FCNC processes seemed to 
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be extremely inhibited in nature (actually forbidden at the tree-level and highly suppressed 
at the one-loop level) the main motivation for the aforementioned GIM mechanism [12], 
nowadays embedded in a natural way into the current formulation of the SM - essentially 
into the unitarity of the CKM matrix. It is remarkable, however, that the degree of 
suppression at one-loop order can vary from one process to another in a dramatic manner. 
For instance, in the b-quark sector the radiative B-meson decay has a branching ratio 
B{h — > S7) ~ 3 X 10~^ which, although small, it has been measured experimentally [13] 
with quite some accuracy and it is used in practice to constrain models of new physics. 
In contrast, the FCNC top quark decay t — > eg becomes radically inhibited in the SM, 
13{t — > eg) ~ 10^^^, namely down to limits far below ever being possibly observed [14-16]. 
Amazingly, the top quark decay into the SM Higgs boson is even more unlikely: B{t — > 
cH) ~ 10"^'' [15, 17]. In all these cases, it is their highly "expected unobservability" 
what provides the natural "signature" for potentially unraveling new physics out of their 
study. In fact, the huge GIM suppression in some rare processes within the SM can be 
significantly softened if one accounts for possible SUSY virtual contributions. For example, 
in the case of the extremely rare top quark decay into the SM Higgs boson one can show 
that if h° is the hghtest CP-even Higgs boson in the MSSM [18], then ^(t ch°) can be 
enhanced 10^'^ times as compared to the SM mode and, thus, bring it to the observable level 
~ 10"'^ [19-^21]. Similar results hold for the Higgs boson decay modes into heavy quarks, 
see e.g. [22-26]. Actually, not only SUSY can help here; other alternative extensions of 
the SM, among them the general Two-Higgs-Doublet Model (2HDM) [18] , predict in some 
cases an enhanced, and often distinctive, FCNC phenomenology [27-34]. Put another 
way: by finding experimental evidence of non-standard FCNC processes we can not only 
enlighten the existence of physics beyond the SM but, in favorable conditions, we can even 
tell the kind of new physics hiding right there ^ . 

In this paper we wish to further explore the FCNC physics of the top quark, but in 
this case we focus on the production of single-top quark final states tc or tc through gluon 
fusion (gg) in pp collisions to take place at the LHC. We denote it by pp(gg) ^ tc + tc. 
While it is true that this process is possible within the strict SM, it proceeds through 
(GIM-supressed) charged-current interactions. We have found the following cross-section 
for this process at one- loop level (see Section Q): 

a(pp(gg) ^ tc + tc)sM = 8.46 x 10"^ pb . (1.1) 

Obviously, it is so tiny that it amounts to less than one event in the entire lifetime of the 
LHC! So it is pretty clear that if this kind of FGNC-generated single top quark signatures 
would ever be detected at the LHC, if only at a level of a few dozen crystal-clear events, 
then the presence of new physics could perhaps be the only valid explanation for them. 
We see that the situation with this production process is very similar to the rare top quark 
decay modes mentioned above; in both cases it is the FCNC physics of the top quark that 
provides the extreme suppression within the SM. However, it should be clear that the top 
quark final states in ( |1.1| ) are a particular class of events within the large variety of single 
top quark processes available in hadron colliders [39-42]. 

^For a review, see e.g. [28] and [35]. See also the interesting flavor mixing studies [36-38]. 
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But this is not the only challenge. The LHC, with all its ability to dig deep beneath 
the physics of the top quark, could perhaps be sensitive to the class of single top quark final 
states associated to FCNC processes, pp(gg) tc + tc, provided of course the underlying 
mechanism could be sufficiently enhanced by some form of new physics capable to boost its 
cross-section up to ~pb level. In this study, we will show that the necessary enhancement 
(which amounts to a factor of roughly lO'^ in the total cross-section) could just come from 
the world of the supersymmetric interactions in the general MSSM. 

Interestingly enough, let us remark that for the FCNC process under consideration, 
pp(gg) — > tc + tc, there is no significant competition between the MSSM and the general 
2HDM because there is no enhancement to speak of from the latter. This is in con- 
trast to the situation with the rare top quark decays mentioned above, where the 2HDM 
contributions are non-negligible as compared to the MSSM ones. Moreover, the direct 
production mechanism pp(gg) tc + tc is substantially more efficient (typically a fac- 
tor of 100) than the production and subsequent FCNC decay of the heavy Higgs bosons 
(A" , H*' — > tc + tc) [32]. In this sense the discovery of a bunch of well-identified tc and/or 
tc events could be strong evidence, not only of new physics, but perhaps of SUSY itself. 
Some of these features were already emphasized in Ref . [43] , where it was presented a first 
self-consistent study of this subject (see also [44]). These references, however, reported on 
the computation of the SUSY-QCD effects only. Other studies can be found in [45, 46] 
under different sets of assumptions. In our case we will continue within the general ap- 
proach initiated in [19], and continued in [43]. It means that the flavor-mixing coefficients 
6ij will be allowed only in the purely left-handed part of the 6x6 sfermion mass matri- 
ces in flavor-chirality space, as it is indeed suggested by standard renormalization group 
(RG) arguments [5,47]. Within this well motivated setup we provide here a full treat- 
ment of the SUSY-EW effects and combine them with the SUSY-QCD ones [43] within 
the general framework of the MSSM. We wish to remark that, in contradistinction to the 
aforesaid studies by other authors, we present our MSSM calculation of pp(gg) — > tc + tc 
in combination with the corresponding MSSM effects on the low-energy b ^ S7 decay and, 
therefore, we extract the single top quark FCNC results only in the region of parameter 
space compatible with the experimental bounds on the radiative B-meson decays. This 
procedure is, in our opinion, a self-consistent approach to the computation of the FCNC 
single top quark signal under study. 

The paper is organized as follows. Section^ is devoted to the general formalism for the 
FCNC processes in the MSSM. In section | we summarize the details of our calculation of 
C7(pp(gg) tc + tc) in this framework. The full numerical analysis is presented in section 
^ leaving section ^ to discuss the results and deliver our conclusions. 

2. Formalism: FCNC interactions in the MSSM 

Apart from the conventional charged-current flavor changing interactions in the SM, the 
FCNC processes in the MSSM are driven by explicit intergenerational mixing terms arising 
from the mass sector of the squarks. For a brief review of this topic, our starting point 
shall be to specify the form of the superpotential, which is the crucial piece of any SUSY 
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theory of particle interactions. In our case we will consider the MSSM with arbitrary soft- 
SUSY-breaking terms. The most general gauge- invariant form of the superpotential can 
be cast in terms of chiral superfields (denoted by a hat) as follows [5-7]: 



Wmssm 



(2.1) 



Indices r, s = 1, 2 refer to the components of the SU{2)l doublets, which are combined in 
a gauge-invariant form through e^s (with ei2 = —£21 = 1). The set of parameters yi, yd 
and yu constitute Yukawa coupling 3x3 matrices in generation space. Although explicit 
generation labels have been suppressed here, they will be introduced at due time. Let us 
notice that in more general SUSY theories there are additional pieces of the superpotential 
inducing violation of baryon or lepton number, but in the MSSM they are set to zero 
because one assumes that the R-parity symmetry holds. 

We also need to settle the piece of the soft SUSY-breaking Lagrangian that takes part in 
the squark mass matrix: 



Csoft = - MqQ*Q - MlU*U - MlD*D 

-.4,. - -1 

+ h.c.. (2.2) 



V2Mw 



cos (3 sin /? 



In this expression, Q stands for the SU{2)l quark doublets, while U,D denote the cor- 
responding singlets. Let us recall that each of the above mass and trilinear coupling 
parameters carries a 3 x 3 matrix structure in the flavor space, although we shall not keep 
track of it explicitly. 

We can now collect the different pieces contributing to the general form of the squark 
mass matrix, which come either from the explicit mass terms in the soft-SUSY-breaking 
Lagrangian ( ^.21 ) or from the couplings triggered by the superpotential (|2.1[) after sponta- 
neous symmetry breaking (SSB) of the EW symmetry. If we arrange all such terms in a 
2-dimensional left-right chirality space, we are left with the following mass matrix: 

2 _ ( ^"^l L + + 2/3(T|^ - sin2 % ) M| m, \ 

ruq Ml^ M? ^ + m2 + cos 2/3 sin^ Ow Ml ) 

(2.3) 

where in the off-diagonal mass terms we have defined MJ^^ = — /i cot (3 and Mf^^ = 
Ad - /itan/3. As usual, tan/3 =< > / < >= V2/V1, with + vl = G~p" j^/l, 
defines the ratio of the vacuum expectation values of the two Higgs doublets giving masses 
to the up and down quarks respectively, while T3 stands for the 3th component of the 
weak isospin of the left-handed quark qi^, and Qq denotes its charge. The non-diagonal 
structure of ( p.3p in the chirality basis requires its diagonalization in order to obtain the 
physical mass-eigenstates in terms of the electroweak (EW) squark eigenstates with 
well-defined SU{2)l x U{1)y quantum numbers. If R^'^^ denotes the matrix rotating the 
qth fiavor, we can diagonalize the mass matrix as follows: R^'^^'^ Ai'^R^'^^ = diag(m?^, m?^). 



Notice that each matrix elements in Eq (2.3) is proportional to the unity matrix Isxs in 
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the flavor space. It is worth realizing, however, that such a trivial flavor structure for 
the mass matrix does not provide the most general realization of the squark mass sector. 
Indeed, in the MSSM we have two fundamental sources of flavor violation. One of them 
just mimics the SM one, namely it consists of the flavor mixing among up- and down-like 
squarks triggered by the charged-current interactions induced by the charged gauge bosons, 
the charged Higgs bosons and the charginos. The second one is qualitatively new and is 
caused by the so-called misalignment between the rotation matrices that diagonalize the 
quark and squark sectors or, in other words, the fact that the squark mass matrices in 
general need not diagonalize with the same matrices as the quark mass matrices [47-50]. 
This is reflected in the existence of the gaugino-fermion-sfermion interactions mediated 
by gluinos (g°; a = 1,2,. ..,8) and neutralinos {Xai o; = 1,2, ..,4). Consider e.g. the 
gluino-quark-squark interactions 

jC^gg = -iV2gsT { ul^ V{u),, {T'')UL^ + V{d)ij {T'')dL, ] + h.c. , (2.4) 

with 

V{u) = B\ul)A{ul) , V{d) = B\dL) . (2.5) 

Here z,j are generation indices, T°- are the SU{2>)c generators, and A{ul)-,B{ui) and 
B{dL) are rotation matrices in generation space which relate the electroweak and the 
mass-eigenstates; e.g. A{ul) rotates up-quarks and B{ui) rotates up-squarks, etc. Notice 
that in the down sector we only need to rotate squarks through B{dL) because after SSB 
of the gauge symmetry the down quark matrix is already diagonal. This follows from the 



fermion mass matrix structure that emerges from the superpotential (2.1) in generation 
space after the Higgs bosons acquire VEV's and spontaneously break the EW symmetry. 
Let us consider only the quark sector. 





-Cm, = -qi, , ^ ] {yd)ij dn^ - qli \ „ {yu)ij urj + h.c. (2.6) 



We can rotate qi- and dji- in generation space until the mass matrix for down quarks, 
ij, becomes diagonal, but then the mass matrix for up quarks, (t>2/\/2) {Uulij, 
will in general be non-diagonal because qi- was already rotated. By inspecting the charged 
current interaction of quarks, this immediately implies that the ordinary CKM matrix is 
just Uqy.m = A'^{ul). Similarly, from the charged current for squarks we read off the 
corresponding CKM matrix in the squark sector: C/scKM = B^{uL)B{dL). Therefore one 
finds a relation between the CKM and SCKM matrices: 

C/SCKM = V{u) UcKM V{d)^ . (2.7) 

As a result, in the MSSM we need three unitary matrices to parametrize the flavor chang- 
ing interactions, one is the ordinary CKM matrix and the other two are associated to the 
new FCNC gaugino-quark-squark couplings (p.4|). Neglecting this second source of flavor 
changing interactions (i.e. assuming that the matrices V{u) and V{d) are unity in flavor 
space) would be tantamount to assume that quarks and squarks diagonalize simultaneously. 
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i.e. [/sCKM = UcKM- This is the super-CKM basis approach to the FCNC processes; it 
assumes that these processes appear at one-loop only through the charged current interac- 
tions (from W^, charged Higgs bosons and charginos and with the same mixing 
matrix elements as in the Standard Model CKM matrix. However, in general we expect 
that the two sources of FCNC should be active in the MSSM and we will take them both 
into account in our calculation. 

These observations turn out to be crucial for the discussion of the flavor-changing 
processes in the MSSM because it means that we can extend the simple 2x2 squark 
mass matrices in chiral space into 6x6 mass matrices in (flavor) (g)(chiral) space. We shall 
comment below on how to parametrize the flavor mixing terms. For the moment we note 
that, due to the aforementioned flavor mixing, the squark mass matrix diagonalization 
process must be extended as follows: 

6 

C = E^i?«f' (« = 1,2,...,6) 
6=1 

i^Wt^^l^W = diag{m|^,...,m|} {q = u, d) , (2.8) 

where A^^. are the 6x6 square mass matrices for squarks in the EW basis {q'^), the eigen- 
values being denoted rn-^. Indices run now over 6-dimensional space vectors with suitable 
identifications. For example, for up-type squaxks a = 1, 2, 3, . . . , 6 = ul, ur, cl, ■ ■ ■ , tn, and 
a similar assignment for down-type squarks. Furthermore, let us notice that the SU{2)l 
gauge invariance of the MSSM Lagrangian imposes certain restrictions over the up-squark 
and down-squark 6x6 soft-SUSY breaking mass matrices, specifically in their LL blocks, 
as follows: 

{Mi)^^ = K{Mi)^^KK (2.9) 

where K stands for the CKM matrix (previously denoted Uckm for convenience). It 
is thus clear, in particular, that both squark matrices cannot be simultaneously diagonal 
(unless they are proportional to the identity) and, therefore, they cannot be simultaneously 
diagonal with the up-like and down-like quark mass matrices either. This is again a reflect 
of the misalignment effect between the mass matrices of quarks and squarks. 

Despite what we have just argued above, within the context of Grand Unified Theories 
(cut's), one usually assumes that the parameters should be aligned at the character- 
istic high energy scale Mx ~ 10"*^^ GeV of these theories (that is to say, the quark and 
squark mass matrices should diagonalize simultaneously at Mx)- But even within such 
theoretically-motivated scenario, it can be shown that the renormalization group running 
of these parameters down to the EW scale would again destroy the primeval aligned con- 
figuration [5,47,48]. It is therefore wiser to take the misalignment into account right 
from the start in the calculation. The most common way to parametrize it is by defining 
the following dimensionless quantities, 5/^^, being A,B = 1,2 the chirality indices and 
i,j = 1,2,3 the flavor ones, in such a way that we can set the non-diagonal squark mass 
matrix elements to be: 

{Mf^^) = Sf^^rhfrhf (z/j), (2-10) 
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where ihf stands for the soft-SUSY breaking parameter of a given chirahty and flavor. (No 
sum over repeated indices here.) It is very common to set all the mass parameters rhf 
equal to a generic SUSY scale MguSY- 

As far as we are dealing with FCNC processes involving the top quark, the most rele- 
vant mixing parameters are those ones relating the heavy up-like flavors among themselves, 
thus essentially t — c transitions parametrized by (523 ("u) above. In close relation to them we 
have the b — s transitions controlled by the parameter {d) ■ Only these mixing parameters 
are expected to be large in GUT's and, moreover, they are not significantly constrained by 
phenomenological considerations. The experimental bounds on the various mixing para- 
meters are derived from the absence of low-energy FCNC processes, which mainly involve 
the first and second generations. For instance, the measurements of the mass splitting in 
K - K° and D - W phenomena [49, 50]. 

Regarding heavy flavors, the phenomenological constraints come from the branching 
ratio of the radiative B-meson decay Bexp(^ §7) and also from the mass splitting in B— 
mixing effects. Clearly, such two processes can only be sensitive to the down-like heavy- 
flavor mixing parameter, 623(d). However, they can also provide information on the allowed 
values for 523{u) since both up and down-like flavor-mixing parameters must necessarily 
be related through the SU{2)i symmetry ( |2.9D . As advertised, in our framework we limit 
ourselves to consider flavor mixings only in the LL-block of the squark mass-matrices, the 
only ones which are well-motivated by RG arguments. Thus, the relevant piece in our 
calculation will be the LL sector of the up-type squark 6x6 mass matrix, which can be 
rewritten in the following manner: 



m)LL = ^SUSY 1 523{u) I . (2.11) 

LL 

Similarly for Mj with 523{u) S23{d). Squark mass-eigenstates follow from diagonalization 




of these matrices through Eq. ( |2.8| ). The mass-eigenstates of (^) are recovered by setting 
the mixing parameters 6 to zero, as could be expected. 

Once we have discussed where the flavor-mixing source is rooted in the MSSM, we 
must now trace back its role at the Lagrangian level [51]. The misalignment between 
the diagonalization matrices in the quark and squark mass sectors triggers the presence 
of couplings of the guise gluino-quark-squark (in the SUSY-QCD part) and neutralino- 
quark-squark (in the SUSY-EW one) that allow the interaction of quarks having the same 
charge but belonging to different generations. At the 1-loop level it is also possible to have 
this kind of flavor-changing interactions mediated by the ordinary SM charged currents, 
but in the i2-odd part of the MSSM we also have the chargino-up-quark-down-squark 
interactions and the charged Higgs-up-quark-down-squark vertices. In all such cases the 
SUSY nature of the couplings allow the resulting process to bypass the SM GIM mechanism 
and provide non-suppressed FCNC events (the charged Higgs piece is an exception, as 
we shall see). The importance of such effect is correlated with the choice of the MSSM 
parameters, in particular those specifying the soft-supersymmetry breaking and, of course, 
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the explicit intergenerational mixing 5^^, which are the most relevant ones for the flavor- 
changing dynamics in the MSSM. For the SUSY-QCD coupling one must work out the 
super symmetrized gauge interaction piece: 

^A^^ = -i^9srk>^%T'')ki^i + h.c, (2.12) 

where T'^ are the gauge group generators, the indices k, I denote the corresponding gauge 
quantum numbers (color, weak isospin) of the interacting particles and A'^ stands for a 
generic gaugino fleld. To extract the FCNC vertices one must include the generation in- 
dices in these interactions. For the particular case of the gluino-mediated interactions, this 
was done in ( |2.4| ). For the practical calculations we will use the extended 6x6 diago- 
nalization matrices defined in ( |2.8D . Therefore, by plugging the squark mass-eigenstates 
in this expression we can rephrase the result in the mass-eigenstate basis and in terms of 
four-component Dirac spinors (for both the gluino and quarks). In the up quark-squark 
sector we get 

6 3 

C^uu = -iV2gsY,Y.<{^aI*tPL-Ri%l-,)gPR)n, + h.c, (2.13) 

a=l 6=1 

and similarly for the down quark-squark sector. Here we have omitted color indices for 
gluinos, quarks and squarks. Notice that while the sum over index a runs over the whole 
flavor (8) chirality space, index b runs only over generations because we are already using 
the standard projectors Pl,r = (1/2) (1 =F 75) to set the chirality of the quarks. 

A similar analysis can be performed to obtain the corresponding Lagrangians describ- 
ing the flavor-changing interactions in the SUSY-EW sector. The calculations are slightly 
more involved since such terms arise from the combination of the SUSY-gauge piece ( 2.12| ) 



together with the higgsino-quark-squark Yukawa couplings dictated by the superpotential 
(cf. Eq. (|2.1| )). Moreover, because of the EW symmetry breaking, the higgsinos and gaug- 
inos mix together to give the final physical eigenstates, the neutralinos Xa i'^ — 2, 4) 
and charginos = 1)2). We shall quote here the final result for such interaction 

Lagrangians (a detailed derivation can be found in [6] and references therein) . For the case 
of the neutralinos, we get: 



4 6 3 



a=l a=l 6=1 



+yuR^:^*N^APR-^ tanewR^:il,^K,Pn 



Ub + h.c, (2.14) 



where g is the weak SU{2)l gauge coupling constant. A few words about notation: in- 
dex a is running over the four neutralino states, Naf3 being the diagonalization matrix 
that provides the neutralino mass-eigenstates, N*Mj^gN^ = diag{m^^^^,m^Q2^''^xo'ii'^xoi)^ 
while Uu stands for the corresponding Yukawa coupling and 6w is the weak mixing angle 
{e = g sin^v^)- Similarly the chargino-up-squark-down-quark interaction Lagrangian can 
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be cast in the following form: 



2 6 3 3 



Xud 



0=1 a=l 6=1 c=l 



+yd R'^' C//32 Pl - Vu RZ' V^, Pr Kl + h.c. 



(2.15) 



This time the standard CKM matrix also needs to be taken into account because of the 
charged-current mixing between up-like squarks with down-like quarks. Again, ?7, V refer to 
the diagonalization matrices of the chargino mass, such that U*M^±V'^ = diag{m^+ ,m^-). 

3. Single top-quark production through FCNC processes in the MSSM: 
computation procedure 

In the following we will concentrate on the analysis of the single top-quark production 
by direct supersymmetric flavor-changing interactions at the LHC, namely the processes 
leading to tc or tc final states. The leading mechanism is the gluon fusion channel: 
pp(gg) — > tc -|- tc (see Section ^ for a full list of Feynman diagrams). It should be clear 
that (T(pp(gg) tc -|- tc) = 2cT(pp(gg) — > tc). It was already shown in [46] that the 
gg partonic channel largely dominates over the qq one at the LHC. Although there are 
previous studies of this process in the literature within the MSSM and adopting different 
approximations [43,45,46], a closer look is highly desirable from our point of view. This is 
so because the kind of simplified assumptions made in some of the previous analyses do not 
shed sufficient light on the possibility that this process could be sufficiently enhanced in 
the MSSM as to be considered realistically at the LHC. We will comment on the differences 
among these approaches later on. In the present paper we carry out our calculation within 
the framework of [43,44], which was first delineated in [19]. 

Throughout the present work we have made use of the standard algebraic and numerical 
packages Feynarts, FormCalc and LoopTools [52-54] for the obtention of the Feynman 
diagrams, the analytical computation and simplification of the scattering amplitudes and 
the numerical evaluation of the cross section (up to the partonic level). Notice, however, 
that we also need to address the computation of the total hadronic cross section in order 
to account for the physical process, a pp collision, to take place at the LHC. To this 
aim we have made use of the program HadCalc [55] ^, while several cross-checks have 
also been performed with other independent codes implemented by us. Throughout our 
calculations we have settled both renormalization and factorization scales at a common 
value, chosen to be half of the production threshold = fip = \{nT'c + rnt). Concerning 
the parton distribution functions (PDF's) involved in the long-distance dynamics of the 
hadronic process, we have included the recent CTEQ6AB data set [56] provided by Les 
Houches Accord Parton Distribution Functions Library (v. 5. 2) [57]. 

The computation of cr(pp(gg) ^ tc -|- tc) in the MSSM is not straightforward. It 
involves a number of subtleties that must be carefully handled. To begin with, we must 

^The source code is available on request from the author. 
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deal with the PDF of a gluon, which exhibits a huge slope in the low-momentum region. To 
that purpose, we have implemented a logarithmic mesh for the integration over the partonic 
variables instead of the linear mesh that is provided by default. We must obviously pay 
the price of adding the corresponding jacobian piece to the original integral. Furthermore, 
a double call to the integration subroutine (viz. the Vegas routine provided by the Cuba 
library [53]) has been implemented. The first call provides an adapted grid for the second 
one, in such a way that the convergence is much faster. As a result a good and reliable 
numerical accuracy is achieved (meaning that the values always remain of order 1). 
Last but not least, a second non-trivial subtlety emerges from the fact that one of the 
final states, the c-quark, has a very small mass when compared to the s/S value of the 
scattering process. We are thus very close to a coUinear-divergence regime. Although there 
is no analytical divergence, the mass of the c-quark is low enough to trigger instabilities in 
the code when integrating over very small angles. We have carefully studied the problem 
and have included a tiny angular cut (chosen to be such that sin^ 9 < 0.03) in order to avoid 
the aforementioned instability. We have checked the dependence of the final calculation 
on the choice of the angular cut (the total hadronic cross section can change about a 15% 
when moving from sin^ 9 < 0.01 to sin^ < 0.1), thus no dramatic changes occur when 
tuning the cut within reasonable ranges. 

In regard to the calculation of the amplitudes contributing to the relevant process 
under consideration, pp(gg) ^ tc -|- tc, let us note that the leading order is the 1-loop 
level. This is a common feature when studying FCNC processes in any renormalizable 
theory (due to the lack of FCNC tree-level interactions). This implies that one need not 
renormalize the bare parameters nor the Green's functions as there are no explicit terms in 
the interaction Lagrangian where to absorb the UV divergences. In other words, the overall 
amplitude of the process should be already finite as soon as we add up all the diagrams 
contributing to that process. In order to check the finiteness of the resulting amplitude, 
we have made use of a standard numerical procedure provided by FormCalc. 

In the following section we present our final numerical results. We shall not provide here 
analytical details of the complicated algebraic structures appearing in the calculation of the 
many one-loop diagrams involved (see Figures 1-5). We have carried out the computation 
in a fully automatic fashion by means of the numerical and algebraic tools mentioned 
above, and of course we have previously submitted our codes to many important tests and 
non-trivial cross-checks of different nature. 

The calculation of cr (pp(gg) tc + tc) has been linked to the one-loop calculation 
of B{h — > S7) in the MSSM, so that after enforcing this low-energy observable to stay 
within the experimental bounds we have obtained the desired cross-section only in this 
allowed region of the MSSM parameter space. Our computation of ;B(b ^ 37) contains 
the complete leading order (one-loop) MSSM computation including the flavor-violating 
couplings. Specifically, we include the contributions to the high energy operators from 
the SM {W^ loops), SUSY-QCD (gluino loops) and SUSY-EW (chargino-neutrahnos and 
Higgs boson loops). The Wilson coefficient expressions have been taken from Ref. [58] 
and they are evolved using the leading order QCD rcnormalization group equations down 
^Reference [58] contains a computation of B{h 37) in the MSSM including some two-loop parts, but 
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to the bottom mass scale. However, at certain stages of our work we use only a part of 
these corrections, this will be clearly indicated in the text below. 

4. Numerical analysis 

To start with, let us present the calculation of the cross-section for the process pp(gg) 
tc + tc within the context of the SM. The Feynman diagrams describing the interaction 
at the partonic level in the t'Hooft-Feynman gauge are shown in Fig. ||. In this gauge, 
the covariant sum over the polarization states of the gauge bosons yields the relation 
e*{k, A) e{k, A) = —g^u- Notice that in the present situation it is unnecessary to intro- 
duce the Faddeev-Popov ghost-field contributions since the current process involves only 
external gluon lines and it thus suffices to restrict the above sum to the two physical degrees 
of freedom carried by the gluons. This is straightforwardly done within the framework of 
the standard computational tools of Ref. [52-54] and allows us to get rid of the spurious 
modes of the quantized gluon field. Also worth emphasizing is the effect of the GIM sup- 
pression, which is associated to the SM diagrams of Fig. |l[ We take for instance a vertex 
correction diagram driven by the exchange of a charged W boson with a pair of quark lines 
closing the loop, and then sum over fiavors. The result is a form factor of the type 

where Kij denote again matrix elements of the standard CKM matrix, and z is a flavor 
index that runs over the down-like quark states d,s,b. We have also included a standard 
numerical suppression factor from the one-loop integral. The additional GIM suppression 
is of dynamical origin within the SM and it stems from the unitarity of the CKM matrix. 
As a result the overall behavior of the form factor amplitude goes like ~ Gprnf, where 
the rrii correspond in this case to down-like quark masses circulating in the loops. The 
cross-section ate = o"(pp(gg) tc) gets suppressed as ate ~ It is thus not surprising 
that we finally get the value early indicated in Eq. ( |1 . 1| ) , which amounts to about 10~^ fb. 
The order of magnitude of the cross-section at the parton level roughly follows from naive 
power counting and educated guess. Using ( [4.1[ ) we get: 

-(gg-tc)~^(^)^G^m^)^ (4.2) 

where we have included factors from the strong (a^ = Qs/Att) and weak {Gp ~ /M'^) 
interactions. For the estimate we include only the bottom quark contribution (with matrix 
elements Vhe — 0.04 and Vtb — 1) as the other terms are suppressed either by very small 
CKM matrix elements or very light quark masses^. At the LHC energies the previous 
estimate provides the cross-section within the ballpark of the exact result after convoluting 

only the one-loop contributions have been used for the present work. 

*This is similar to the kind of ansatz made by Gaillard and Lee to predict the charm quark mass [59], 
except that here it is the bottom quark that gives the dominant effect because the external quarks are 
up-like. 
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Figure 1: Feynman diagrams corresponding to tlie partonic process gg tc in the SM. Tliere is 
a similar set of diagrams leading to the tc final states; this is implicitly understood here, and also 
in the next four figures containing Feynman diagrams. 



with the parton distribution functions. However, in order to understand the dynamical 
mechanism of enhancement of the SUSY interactions (see below) it will suffice to compare 
with the partonic contribution ( |4.2| ). 



Using the exact (numerically computed) result (1.1), we find that this cross-section 
is literally invisible; even assuming a total integrated luminosity of 1000 fb"^ it amounts 
to one tenth of event during the whole lifetime of the LHC. This result supports quite 
convincingly the idea that the eventual detection of such kind of FCNC processes could 
give us important hints of some form of physics beyond the SM. 



Regarding the choice of SM parameters, we have taken the heavy-quark masses and 
coupling constants given by their corresponding renormalization group running values at 
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Figure 2: Feynman diagrams involving gluino and squark loops for the partonic process gg — *■ to 
in the MSSM. They constitute the SUSY-QCD contribution to the single top quark production in 
the MSSM. 



mcifJ-R) (GeV) 


rubinB) (GeV) 


mtifiR) (GeV) 






sin^ 9w{p-r) 


0.877 


3.024 


183.365 


0.1170 


1/128.89 


0.23 



Table 1: Values of the SM parameters at the scale (see the text). 



the renormalization scale of the process, see Section ^. The running masses and coupling 
constants have been explicitly computed by means of the (5 and 7 functions at the one-loop 
level (in the case of a^, we have made use of an specific subroutine from the CERNLIB). 
The obtained values are displayed in Table 0. 
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Xi 

Figure 3: Feynman diagrams involving chargino and squark loops for the partonic process gg — > tc 
in the MSSM. This subset of diagrams is a part of the charged-current SUSY-EW contribution. 




Figure 4: Set of Feynman diagrams involving neutralino and squark loops for the partonic process 
gg — > tc in the MSSM. It defines the neutral-current SUSY-EW contribution. 

Next we evaluate the SUSY-QCD contribution to o'(pp(gg) tc-|-tc). The optimized 
set of values that we have found for the MSSM parameters is indicated in Table ^. Below 
we give some details on their determination, which basically follows the method of [43]. 
However, here we have performed the calculation with a slightly different set of parameters; 
in particular, the SM parameters have been improved by using their RG running values. 
The SUSY-QCD corrections are driven by all possible 1-loop diagrams (vertex corrections. 
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Figure 5: Feynman diagrams involving charged Higgs loops for the partonic process gg tc. 
Together with the chargino diagrams in Fig. ^ they represent the full non-SM charged-current 
contribution to single top quark production in the MSSM. 



tan P 


5 


AtiGeY) 


2238 


A(GeV) 


2000 


mg(GeV) 


200 


MsusriGeY) 


746 




400 




0.7 



Table 2: Set I of MSSM parameters, that optimizes the SUSY-QCD contribution in the absence 
of SUSY-EW effects. 

self-energy insertions and box diagrams) involving gluinos and squarks (cf. Fig. §). In the 
following we describe the behavior of the SUSY-QCD contribution to the total hadronic 
cross section ats as a function of a given parameter at fixed values of the others, taking 
Table | as a reference. The corresponding results are reported in Figs. 0, 

To begin with, we consider the curve ate as a function of tan p. What we find is that 
the cross section grows steadily until reaching a saturation regime at values of tan/3 ~ 15. 
The shaded region is ruled out by the experimental determination of the branching ratio 
SexpO^ sj) [60]. The excluded region reflects that the MSSM calculation of crtc in it 
yields a value of B{h S7) out of the experimental band allowed for this observable in 
the range i3exp(b — > 37) = [2 — 4.5] x 10~^ at the Su-level ~ see [60] for details. It can be 
proven that the overall MSSM amplitude for b — > S7, and the purely-SM one, both must 
have the same sign [61]. We have included this restriction also in our numerical codes, so 
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Figure 6: SUSY-QCD contribution to tlie total cross section ats (in pb) and the corresponding 
number of events per 100 fb"^of integrated luminosity at the LHC, as a function of a) tan /3 and b) 
At for the parameters quoted in Table ^ The shaded region in a) is excluded by Bexp{i> S7). 



that we automatically enforce the different scannings over the MSSM parameter space to 
be consistent with both the experimental band and the sign criterion. 

We should clarify at this point that up to now we are just retaining the SUSY-QCD 
(gluino-mediated) contributions for the computation of the two observables B{h — > S7) and 
a"(pp(gg) tc + tc). In other words, at the moment we neglect the EW effects both from 
supersymmetric particles and Higgs bosons in all these processes as if they were exactly 
decoupled. At due time we will switch them on in combination with the SUSY-QCD effects 
to evaluate the full MSSM result. 

It is worth noting the strong dependence of cite on the trilinear coupling At (cf. Fig. 
|6|b), where we have included the approximate constraint \At\ < 3MsusY to avoid color- 
breaking minima. We see that ate changes around two orders of magnitude along the 
explored range. The dependence of the cross section as a function of the SUSY-breaking 
scale (Msusy) and the gluino mass (rrig) is given in Figs. ^,b. In both cases ate decreases 
with the mass scale, as expected from the decoupling theorem, but this feature is more 
accentuated with the parameter rrig. For instance, ate becomes 10 times smaller when 
increasing the gluino mass from rrig = 200 GeV to rrig = 500 GeV. For the given values of 
the parameters in Table |2|, MsijSY cannot be smaller that the value indicated there, the 
reason being that we must respect the lower mass limits on the squark masses. For the 
latter we just take the LEP limits. This means that we do not exclude from our scanning 
regions of the parameter space where some physical squark masses can be as light as 90 
GeV, hence we assume rriq > 90 GeV in our analysis^. 

There is also a monotonously decreasing trend when scanning over the higgsino mixing 
parameter fi (Fig. ^) , although in this case the variation involves less than one order of 

^For a general overview of the different strategies and up-to-date results concerning the squark mass 
bounds, see Ref. [13] 
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Figure 7: SUSY-QCD contribution to tlie total cross section ats (in pb) and the corresponding 
number of events per 100 fb"^of integrated luminosity at the LHC, as a function of a) Msusy and 
b) rrig for the parameters quoted in Table |^. 

magnitude in the allowed range. Concerning phenomenological bounds, values of > 
800 GeV are excluded by the observable Bexpi^ — *■ S7), and also because of the bounds on 
the lightest squark mass. LEP bounds also exclude < 200 GeV since otherwise the 
chargino mass limit, ?7z^± > 94 GeV, would be violated. 

The most dramatic dependence of the SUSY-QCD contribution to ate arises from the 
explicit flavor-mixing terms. This can be seen at work in Fig. ^3, where we scan over S2^{u). 
In this case the cross section grows from (the gluino-quark-squark coupling vanishes in the 
limit (u) = 0) to 2(7tc ~ 1 pb at the maximum allowed value of the flavor-mixing, viz. 
^23 i'^) — (s^^ below). Larger values of (u) are excluded by the lower experimental 
limits on the squark masses (we recall that the flavor-mixing parameters participate in the 
diagonalization of the squark mass matrix, see Eq. ( p. 111 )). 

It is important to realize that 62^ {u) is also constrained by BexpO^ — ^ S7). This was 
advanced in Section |2|. As we shall next argue, the Set I of MSSM parameters (cf. Table |2[) 
does maximize the SUSY-QCD contribution to ate within the present phenomenological 
constraints on the MSSM parameter space. We will refrain from writing cumbersome 
analytical expressions for the exact formulas. However, we can provide the main analytical 
ingredients of the calculation in a schematic way as they will be useful to understand the 
physical origin of the SUSY enhancements. Let us illustrate the procedure by calculating 
the approximate optimal value of 62^ (u). The starting point in this discussion is the general 
form of the SUSY-QCD contributions to the cross section. It will suffice to consider the 
partonic cross-section since all the distinctive dynamical features are already contained in 
it. From the formulae of Section^, and educated guess, we find 

^ . M!€ ^?(.4.- M/tanfff 

Let us briefly explain the origin of some terms in this expression. There is a loop numer- 
ical factor as in ( [4. 2D . All interaction vertices are strong and therefore we must have a 
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Ug dependence. A bit more subtle is the factor of mt {At — fi/ tan/3) at the level of the 
amplitude, which originates from the chirality flip of the gluino in the loops and the cor- 
responding chirality LR transition of the top-squark. This produces a factor of Mj^j^ (in 
the amplitude) emerging from the stop mass matrix - see Eq. (O) - which goes square 
in the cross-section. There is of course also the (square of the) important flavor mixing 
factor 623{u) stemming from ( p.ll| ). Finally, the factor Mg^gy the denominator of the 
amplitude (appearing as Mg^gy in the cross-section) is related to the SUSY particles cir- 
culating in the loops. As a matter of fact, this factor may get contributions from squark 
masses, the gluino mass or some combination of them. By comparison of equations ( [4.3| ) 
and (|4.2| ) we can easily deduce the expected order of magnitude of the SUSY enhancement 
factor in the total cross-section. For simplicity, let us approximate M|^^ ~ mtAt, which 
satisfies the conditions of Set I of parameters in Table ^. Then the ratio between the SM 
and SUSY-QCD partonic cross sections yields 



a{gg tc) 



SM 



I He 



SUSY-QCD 



mt 



mt 



mt 



(4.4) 



Using the inputs on Table § we can easily check that the SUSY-QCD cross-section is of 
order 10'' — 10^ times larger than the SM one. This result will essentially persist after 
convoluting with the parton distribution functions. Therefore, the cross-section will move 
some 7 orders of magnitude upwards, namely from the tiny value (^]^) up to values of order 
1 pb, which is perhaps sufficiently large to have a chance of being detected. 

In a similar way, we can easily determine the general form of the branching ratio 
B{h — > S7) in the MSSM. In this case, we will only indicate the SUSY factors relevant for 
our considerations. To start with, a chirality flip in the bottom-squark line is involved. 
This follows from the structure of the Wilson operator describing the effective low-energy 
interaction [58]. Thus, a factor (A^ — /itan/3)'^ arises in the corresponding branching ratio. 
The result is the following: 



B{h S7) ~ \6., 



231 



2 ml{Ah — fitan j3y 
^•^SUSY 



(4.5) 



where the parameter (523(6) emerges from a mass matrix similar to ( p.llj ) but for the 
bottom-squark sector. In this respect, let us recall the relationship (|2.9| ) linking the LL 
blocks of the stop and sbottom mass matrices, the only ones involved in this calculation. 

mt{At — ^/tan (3) / M'g^gYi which is nothing but the result of 



Let us define S^^^ 



rewriting the off-diagonal LR matrix element of Eq. (2.3) in terms of the flavor- mixing 



parameters defined in Eq. ( 2.11 ). With this definition, we can think of the curves having 
a constant value of ate as hyperbolae in the ^33^ — 
squark mass matrix in the following simplified form: 



plane. Next we take the up-type 
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(4.6) 
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Figure 8: SUSY-QCD contribution to tlie total cross section ats (in pb) and the corresponding 
number of events per 100 fb"^of integrated luminosity at the LHC, as a function of a) ^ and b) 
523^ (u) for the parameters quoted in Table ||. The dashed area in b) is ruled out by the lightest 
squark mass bound. 



Upon diagonalization, we can see that the allowed values of the squark masses are such 
that they fullfill the equation {623)'^ + ((^la^)^ = (1 - "^|i/^|c/5y)^i which is nothing but a 
circle of radius R = 1 — "i^^-^ / ^su SY ■ Taking advantage of the above geometrical picture, 
the only thing we have to do in order to determine the maximum of ate is to identify the 
particular point of the straight line = §2^, in the — plane, for which the 
outermost hyperbola ate =const. is tangent to the circle of radius R. Thus, we find that 
the (approximate) value of 62^ that maximizes ate reads: 



^23 



V2 



1 + 



1 + 



2 

2 91 

9 nit'^ 



1/2 



0.7, 



(4.7) 



where has been taken 90 GeV. For comparison, 62^ ~ 0.68 if niq-^ = 150 GeV, so it 
does not critically depend on the lower bound on the squark masses. The parameter value 
(4.7) is precisely the one quoted in Table ||. The corresponding value for the total cross 
section is o"(pp(gg) — > tc + tc) ~ 1 pb, which translates into a production rate of around 
10^ events per 100 fb"^of integrated luminosity. Notice that the Set I of MSSM parameters 
is a convenient choice to maximize the cross-section since it satisfies A}, — fi tan /? = and, 
as a result, the constraints imposed by BexpO^ — *■ §7) are automatically satisfied and impose 
no additional restrictions on Si2'i{d) and, hence, not on Si2'i{u) either. 

Besides the b — > S7 bounds, it is also of interest to explore the impact of the recent 
double-side constraints on the mass splitting between the B° and states [62], which 
yield (at the 90% CL): 

ITps'^ < A Ms < 21 ps'^ (4.8) 
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Figure 9: — constraints in the 623" {u) — rrig plane. The present computation takes into 
account the SUSY-QCD contributions to the B° — b" amphtude within the mass insertion appro- 
ximation (See e.g. [49]. ) 





m-g (GeV) 


a{ph) 


0.70 


791 


0.00722 


0.6 


731 


0.00422 


0.50 


656 


0.00332 


0.40 


556 


0.00324 


0.30 


425 


0.00388 


0.20 


236 


0.00732 


0.18 


200 


0.00806 



Table 3: Total cross section cr (pp ^ tc + tc) (in pb) for the lowest allowed values of the gluino 
mass rUg and corresponding [u) values in Fig.|| obtained from the B° — B^ bounds. The rest 
of MSSM parameters are fixed as in Table ||. 



The former result can be rephrased as follows (see [63] for details): 

MsusY 



0.55 < 



1 + 



M 



SM 



< 1.37, (4.9) 



where Ai = (B^I'H^'^^^ |B^) stands for the transition matrix element that describes the 
— Bj, mixing, in which supersymmetric radiative corrections may certainly play a role. 
A number of studies [64-66] have pointed out the possibility that those diagrams involv- 
ing gluino exchange could carry the bulk of the SUSY contribution to A^. It has also 
been shown that the aforesaid SUSY corrections are very sensitive to combinations of 
flavor-mixing parameters of different chiralities, namely S2^{d) '^23^(d), ^23^{d) 62^{d). Let 
us recall, however, that we are only allowing non- vanishing LL-flavor-mixing parameters, 
^23 7^ 0- Therefore, for consistency, we will study the influence of the B^ — B^ mixing 
bounds under this same assumption. The effective AB = 2 Hamiltonian describing the 
B^ — B^ mixing can be expressed in terms of local operators and the corresponding Wilson 
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coefficients, see e.g. [64] for a result obtained within tlie so-called mass insertion appro- 
ximation (in which the sfermion propagators are expanded in powers of the flavor-mixing 
parameters up to the linear order). Focusing on the SUSY-QCD (viz. gluino-mediated) 
contribution, and setting to zero all the flavor-mixing parameters but (^23^ (u) (and thus 
(5|3^(d) through the SU{2)l relation (|2j)), we are left with a single Wilson coefficient for 
the MSSM contribution, 



,2 



^1^/^) = -^^ 24x/6(x) + 66/6(x) {6^i^id)f. (4.10) 
216m| L J 



This coefficient must be computed at the energy scale where the actual mixing takes place 
(/i = TJih) by means of the corresponding RG equation (Cf. Refs. [49,64]). Finally, we may 
compute the amplitude for the — mixing process, 

/ttO, 



1 

~3 



lcfinib)mB,JlBiip), (4.11) 



where we have used 



(B^IOiIbO) = -^niBjlBiii,), (4.12) 

Bi being a "bag factor" (expected to be of order ~ 1), whose numerical value is determined 
by lattice techniques [67]. 

A glimpse at equations ( |4.10|) and ( 4.11 ) suggests that the B^ — B^ mixing bounds 



may impose further restrictions on rrig and 62^. However, we also notice that a "bag" 
dependence enters the analysis of this system. Namely, the coefficient Bi is a clear signal 
of strong binding effects that cannot be accounted for in the perturbative framework. 
The appearance of these coefficients is characteristic of calculations involving bound state 
systems by strong interactions. One has to resort to lattice calculations to estimate them. 
In this sense, this kind of calculations are imbued of higher theoretical uncertainties. Since 
these bag parameters are not present in the calculation of ;B(b — > S7), we deem more suited 
to treat the respective constraints derived from b ^ S7 and from B^ — B^ at different 
levels of applicability, the former being more solid than the latter. This way of proceeding 
also helps to more clearly differentiate the impact of the different kind of constraints on 
our FCNC production process. 

Plugging the above expression into our codes, we can single out which regions in the 
nig — (u) plane are excluded (or allowed) by the experimental bounds of Ref. [62]. The 
resulting plot is presented in Fig. ^. As could be expected, only heavy gluino masses can 
now be compatible with large flavor-mixing values, whereas light gluino scenarios enforce 
(523^ (n) to stay low. The maximum cross sections of our FCNC process for different values 
of TUg and 523" are presented in Table |3[ Notice that the optimal regimes are attained at 
the edges of such range, that is to say, when 62^ (u) is maximum and rrig is as light as 
permitted by the experimental bounds on the B^ — B^ mixing - analogously, when rrig 
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is minimum and (u) is as large as allowed. In such optimal regimes we are left with 
maximum cross sections barely of the order of <t ~ 10"2 pb - implying ~ 1000 tc(tc) 
pairs per 100 fb"^ of integrated luminosity. Interestingly enough, such non-negligible rates 
could perhaps be attained even for very heavy gluinos {rrig ~ 800 GeV) in the region of the 
highest allowed value of 62^, given by (4.7). However, we should stress that, in this range 
of values of the flavor-mixing parameters, the validity of the mass insertion approximation 
could fail significantly. This means that, in order to estimate with some reliability the 
allowed/excluded region of Fig. say for 62^ > 0.5, it would require the inclusion of 
additional terms in the expansion of the sfermion propagators, or simply to go beyond that 
approximation through an exact calculation of the flavor- mixing effects at one-loop. This, 
more refined, calculation is at the moment not available in the literature and it certainly 
goes beyond the scope of our main aim in this work. On the other hand, as we have seen, 
there are other sources of uncertainty that could challenge the real usefulness of performing 
such calculation. 

We conclude that the inclusion of the — mixing constraints could further limit 
the SUSY-QCD enhancement capabilities by roughly 2 orders of magnitude with respect 
to the optimal regime in Table |2|, but only if one can apply the restrictions obtained from 
the mass insertion approximation for the range of 62^ near the optimal value ( [4.7D . Recall 
that the optimal regime (where the relevant FCNC cross section can be of order of Ipb) 
was obtained in agreement with the experimental restrictions imposed by b ^ S7 alone. 
However, as warned above, the b — > S7 and B^ — B^ constraints are of different nature 
and, in our opinion, they should not be placed on equal footing. The main reason for this is 
that, for the inclusive decay B Xg^, the theoretical prediction for the branching ratio can 
safely proceed by using the Operator Product Expansion (OPE) within the Heavy Quark 
Effective Theory. Uncertainties related to the non-perturbative effects in this case are 
much lower than when we are forced to follow other strategies, such as the introduction of 
bag factors - see Eq. ( 4.12| ). The latter must be estimated through lattice calculations and 



suffer in general from a larger amount of uncertainty. This intrinsec source of uncertainty 
is part of the reason why we do not consider necessary for the moment to take into account 
a more detailed calculation of the impact of the B^ — B^ mixing constraints on our main 
process beyond the mass insertion approximation. 

We turn now our attention to the SUSY-EW contributions to the process pp(gg) — > 
tc + tc, which we will also treat in combination with the corresponding effects on the decay 
b — > S7. As we have already discussed in Section ^, the SUSY-EW effects are contained 
in the loop diagrams involving charginos (Fig. neutralinos (Fig. ^) and charged Higgs 
bosons (Fig. |5|). It is again of interest to analize the behavior of these corrections to the 
total cross section as a function of the different MSSM parameters. Although we could 
try a similar discussion pattern as that employed for the SUSY-QCD effects, a qualitative 
treatment of the SUSY-EW effects is far more difficult and cannot be easily summarized in 
terms of approximate analytical formulae as we did for the SUSY-QCD effects. Therefore, 
hereafter we limit ourselves to the corresponding numerical analysis. 

In what follows we adopt a second choice of MSSM parameters. Set II (see Table 
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Table 4: Set II of MSSM parameters optimizing the SUSY-EW contribution in the presence of 
SUSY-QCD effects. 



This choice optimizes the SUSY-EW contribution. The reason should be clear: Set II is 
characterized by a large value of the gluino mass, viz. sufficiently large {rrig > 2 TeV) 
that the SUSY-QCD effects are virtually decoupled. At the same time, the values of the 
electroweak parameters /i, Mi and M2 in Set II insure that some of the charginos and 
neutralinos can be relatively light (within their experimental bounds). The corresponding 



results are presented in Figs. [1^, 11, Throughout these pictures, shaded regions are 
excluded by the constraints imposed by i3exp(b sj), while the dashed areas violate any of 
the current mass bounds on the SUSY particles. The respective LEP limits on the lightest 
neutralino, chargino and squark read: m^o > 46 GeV, ^ 94 GeV and rn-g^ > 90 

GeV. We must also insure that the mass bound on the lightest Higgs boson is respected, 
rriifl > 90 GeV [13]. 

We begin our study of the SUSY-EW effects by considering a scan over tan /3. The 
contribution from the charginos (dotted line in Fig. [To|a) shows a mild growing trend. In the 
case of the charged Higgs contribution, the evolution with this parameter (dashed-dotted 
line in Fig. |l^b) exhibits a bigger slope, but the actual size of the contribution remains very 
small; in fact, not larger than 10~^ pb within the allowed range. The neutralino part slowly 
decreases with tan /5 since the neutralino-quark-squark coupling involves the combination 
CyOqg ~ —\At\ + iJ,/tan(3, for At < 0. The leading SUSY-EW contribution clearly comes 
from charginos. Neutralinos give sizeable, but smaller, effects and Higgs bosons give an 
entirely negligible yield, typically three orders of magnitude below that of charginos (at 
least for the sets of parameters considered in this work). In this sense, what we usually 
call the SUSY-EW contribution (defined as the sum of all the chargino, neutralino and 
charged Higgs boson effects) can be considered virtually identical to the chargino plus 
neutralino contributions. In what follows we will avoid plotting the individual charged 
Higgs contribution. 

Consider next the sensitivity to the trilinear coupling At (Fig. [To|c). From the interac- 
tion Lagrangians ( p. 14 ) and ( |2.15| ), it follows that the neutralino piece is the only one that 
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Figure 10: a) Individual contributions from charginos and neutralinos, as well as the total su- 
persymmetric electroweak effect (indicated by SUSY-EW), to the total cross section ate (in pb), 
and the corresponding number of events per 100 fb"^ of integrated luminosity at the LHC, as a 
function of tan/3, b) Similar as in (a), but for the charged Higgs effects only; c) Individual and 
total SUSY-EW effects as a function of At . In all cases the fixed parameters are as in Table ^. The 
shaded regions are excluded by Bexpi^ — > S7) and the dashed regions are ruled out by the mass 
bounds on the lightest supersymmetric particles. 



is sensitive to At. Fig. [To| c shows that its dependence is quaHtatively similar to that of 
the gluino (cf. Fig. ^b), with a minimum at low At and the largest contributions achieved 
at the highest possible values of \At\. Owing to the b — > S7 restrictions, the relative sign 
between At and the higgsino mass parameter fj, becomes of crucial importance. The favored 
range corresponds to At fj, < 0. For the set of parameters collected in Table it can be 
checked that the leading SUSY contribution to B{h — > 37) is driven by the chargino piece. 
This piece is strongly dependent on fj.. In particular, the chargino amplitude accounting 
for b ^ S7 changes its sign when we reverse fi — > —fi. In our analysis we assume, for 
definiteness, /x > and At < (motivated in part by some preference for the sign fi > 
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Figure 11: SUSY-EW contribution to tlic total cross section ate (in pb) and tlie corresponding 
number of events per 100 fb"^ of integrated luminosity at the LHC, as a function of a) /i and b) 
MsuSY for the parameters quoted in Table ^. Owing to the B{h S7) constraint, the allowed 
range for Msusy in (b) is very narrow (non-shaded area). 



from g — 2 oi the muon ) and scan over the mass parameter fi. In Fig.[ll|a we can see that 
there is a mild increasing behavior of the cross-section for both the chargino and neutralino 
contributions. We have also checked that the alternative choice of signs (/u < 0, > 0) 
does not alter dramatically the leading SUSY-EW effects. 

Next we analyze the behavior of dtc under the variation of Msusy, i-e. the common 
squark soft-SUSY breaking mass scale. Both the chargino and neutralino terms exhibit a 
monotonously decreasing behavior as functions of Msusy (Fig. pT|b). The neutralinos are 
the most sensitive ones and decrease around 3 orders of magnitude until they reach an 
almost saturation regime at a value of Msusy ~ 700 GeV. Such saturation arises from the 
particular contributions of the gaugino-quark-squark couplings in this region of the MSSM 
parameter space. Although not shown in the figure, the cross-section eventually vanishes 
at very large values of Msusy, as required by the decoupling theorem. Actually, a large 
portion of the MsuSY range in Fig. 0=" (viz- Msusy > 288 GeF) lies in the shaded 
region and, hence, is ruled out by the B{b 37) constraint. Once again we see that the 
inclusion of this low-energy observable plays an important role to limit the scope of the 
allowed parameter space. 

Finally, in Fig.|l^,b we explore the response of the SUSY-EW effects to a variation 
with the bino (Mi) and wino (M2) mass terms in the gaugino-higgsino mass matrices. While 
the chargino piece remains constant, because it is unrelated to Mi, both the chargino and 
neutralino contributions decrease significantly with M2. For example, in the case of the 
charginos, and for the range explored in Fig.|l^, it involves a change of more than 3 orders of 
magnitude. Let us note that the reason for accommodating relatively light SUSY particles 



®For the SUSY effects on — 2, see e.g. the excellent review [68] and references therein. Let us clarify 
that the MSSM contribution to this observable depends on the value of some slepton masses which, however, 
play a marginal role in our calculation. 
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Figure 12: SUSY-EW contribution to the total cross section ate (in pb) and the corresponding 
number of events per 100 fb'^of integrated luminosity at the LHC, as a function of a) Mi and b) 
M2 for the parameters quoted in Table ^ 



in the spectrum (charginos, neutralinos and squarks) is because we have to balance the 
important effects from the relatively light charged Higgs bosons (M/^± > 170 GeV, see 
Table ^ in the SUSY computation of ^(b —>■ 37). This explains why only a limited upper 
range of values for MsusY is allowed in the current scenario. 

Our numerical analysis has clearly identified the region of the parameter space where 
the SUSY-EW contributions are most favored and can compete with the SUSY-QCD ef- 
fects. To wit: it corresponds to having very heavy gluinos (of a few TeV) together with 
light neutralinos, charginos and squarks (in particular, stops) within the experimentally 
allowed lower mass limits. This is indeed the motivation for having introduced a second 
choice of MSSM parameters (Set II in Table |^. Still, it may be useful to consider a third 
choice, as we comment below. 

In the first part of our study we have considered the SUSY-QCD effects alone. Af- 
terwards, we have switched on the SUSY-EW effects in a situation where the SUSY-QCD 
contribution is subdominant (viz. characterized by very heavy gluinos). Before closing our 
investigation, it may be interesting to consider a third scenario, namely one in which the 
SUSY-QCD and SUSY-EW effects are both present but the former are dominant. (No- 
tice that even when both of them are favored, the SUSY-QCD effects are usually more 
important.) Following our standard procedure, we have numerically determined the cor- 
responding region of the MSSM parameter space under the B{h S7) constraint. The 
result is collected in the form of the parameter Set III in Table ^. 

In Fig. 13a we plot the cross-section ate for Set III, i.e. the SUSY-QCD favored one, 
as a function of S2^{u). Similarly, in Fig. [l^ we show the yield from the SUSY-EW 
favored case, represented by Set II of Table ^ In the two sets we have the concurrence 
of SUSY-QCD and SUSY-EW contributions, but while the former is characterized by 
relatively light gluinos, charginos and neutralinos, the latter contains light charginos and 
neutralinos, together with very heavy gluinos. 
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Table 5: Set III of MSSM parameters optimizing the SUSY-QCD part in the presence of the 
SUSY-EW contribution and preserving B{h — > S7) within the experimental bounds. 
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Figure 13: Combined SUSY-QCD and SUSY-EW contributions to the total cross section ate (pb) 
as a function of (5^3^ (u) for the choices of parameters that optimize a) the SUSY-QCD term (cf. 
Table ||) and b) the SUSY-EW term (cf. Table |). RecaU that cr(pp(gg) tc tc) = 2atc- 



In both panels of Fig. we display the individual contributions from charginos (dotted 
lines), neutralinos (dashed-dotted lines) and the gluino (full lines); dashed lines denote the 
full SUSY-EW contribution. Finally, the double dotted-dashed line in this figure accounts 
for the total sum of all the 1-loop diagrams. The fact that the full SUSY-EW contribution 
is non-vanishing in the limit ^23 — > is due to the tiny charged Higgs boson piece. In 



Fig. 13 a the SUSY-QCD effects are overwhelming as compared to the SUSY-EW ones; 



the gluino curve is dominant by at least 4 orders of magnitude over the chargino effects 
(the leading SUSY-EW ones). For this reason the interference amplitudes between these 
two kind of SUSY contributions can be neglected in this case. The highest production 
cross-section (namely, the maximum of 2atc) in this case can reach ~ 1 pb, or equivalently, 
it implies ~ 10^ events per 100 fb"^ of integrated luminosity. 

The situation changes dramatically when moving from Fig. |l3|a to Fig. 13 b. The 
latter shows the results corresponding to the Set II of parameters (the SUSY-EW favoring 
ones). The most noticeable feature is that, for intermediate values of the mixing pa- 
rameter (say for {u) > 0.2), the contributions from the various sources of SUSY-EW 
contributions become comparable and are, within this parameter setup, larger than the 
SUSY-QCD effects. Owing to the fact that the chargino interactions involve the charged- 
current chargino-quark-squark coupling, the dependence that we observe on (u) is not 
clean, i.e. it is not a direct one; rather, it is obtained via Eq. ( |2.9| ), which relates the para- 
meters (523^ (u) and 523^ (d) through CKM rotation. We have worked out such relations and 
implemented them in our codes. Since in this second scenario the SUSY-EW effects can 
be comparable to those from SUSY-QCD, it follows that the interference terms between 
both sets of SUSY diagrams cannot be neglected any more. This becomes transparent by 
considering the curve in Figure. 13 b describing the behavior of the overall 1-loop contribu- 
tion. Even for low values of the mixing parameter, the SUSY-EW corrections turn out to 
be significant (for instance, for 62^ {u) ~ 0.4 the full 1-loop curve is enhanced by a factor of 
2 with respect to the SUSY-QCD piece). For 5^^{u) ~ 0.6 and above, the SUSY-EW part 
exceeds by roughly one order of magnitude the SUSY-QCD contribution and becomes the 
dominant effect on ate- Obviously, for heavier gluinos {rrig of a few TeV), the SUSY-EW 
part would be more important and it could remain as the only SUSY source of tc and tc 
events. In this case, the maximum value of 2atc remains still be sizeable as it would amount 
to ~ 1000 events per 100 fb'^of integrated luminosity. 

Some observations on previous work on this subject are now in order. The kind of 
systematic analysis presented here was not done in previous work on this subject [45,46]. 
Only some particular regions of the MSSM parameter space were singled out in these 
references and, moreover, the important restrictions due to BexpO^ §7) have not been 
taken into account. Thus e.g. in [46] the bulk of the contribution was missed and the 
authors tended to emphasize that the main effects stem from the mixing in the LR sector 
of the full mass matrix, A^|, namely from 7^ j). Similarly, in [45] the various chiral 

mixing effects are taken into account and, again, emphasis is made on scenarios where the 
maximal effects are obtained from LR mixing. However, it is not obvious to us whether the 
important SU{2)l relation ( p.9D was enforced in this reference when dealing with the LL 
sector. Furthermore, the simultaneous analysis of the Bexp{b S7) within the MSSM is 
also absent in [45] . We have nevertheless tried to compare our numerical results with those 
presented in the latter reference by adapting our codes to accommodate specific features of 
their calculation. Basically, it involves a different choice of SM and MSSM parameters, a 
different criterion for the renormalization scale and the inclusion of some phase space cuts. 
When doing so, we are able to successfully reproduce the results provided in Table II of [45] . 
We should also mention that in [36] an analysis is made of FCNC decay modes and the 
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single top quark production process in which the b ^ S7 constraint is considered. However, 
our main emphasis as weh as the conditions under which our results have been obtained, 
are different. To put it in a nutshell, in contradistinction to all previous approaches, what 
we have shown here goes well along the lines of Ref. [43, 44] (and fully generalizes their 
main conclusion for the complete set of MSSM quantum effects), to wit: in order to get 
sizeable SUSY contributions to cr(pp(gg) tc + tc) (i.e. cross-sections of order of ~ Ipb) 
compatible with the low-energy FCNC constraints, it suffices to retain the LL box of the full 
flavor-mixing matrix, which is, on the other hand, the only one well motivated by standard 
RG arguments. If, in addition, the B° — B° mixing effects are taken into account in this 
same kind of scenario, the maximum cross-section could be further reduced. However, for 
small values of 62^, in which the mass insertion approximation can be relied, the reduction 
is not very significant and applies to a value of the cross-section which is already small 
(~ O.Olpb). 

5. Discussion and Conclusions 

In this work we have studied the single top-quark production by strong and electroweak su- 
persymmetric flavor-changing interactions at the LHC. We have concentrated on the leading 
gluon-gluon fusion mechanism, pp(gg) ^ tc -|- tc, and restricted within the context of the 
Minimal Supersymmetric Standard Model. Among the super symmetric flavor-changing 
interactions in the MSSM, we have the charged current ones induced by charginos and 
charged Higgs bosons, but also the Flavor Changing Neutral Currents (FCNC) triggered 
by gluinos and neutralinos. A self-consistent approach to this high-energy calculation can 
only be achieved by taking into account the important experimental constraints on low- 
energy FCNC processes, essentially on the radiative B-mcson decays (b S7). To our 
knowledge, the first study of the supersymmetric single top quark production processes at 
the LHC under these conditions was performed in [43,44]. Other studies of single top quark 
production in the literature [45,46] simply obviated the implications on that low-energy 
process, which is experimentally well measured and highly sensitive to the dynamics of the 
SUSY interactions [58,61]. However, the study of [43] was confined to the computation of 
the SUSY-QCD part. Here we have extended it by computing the full SUSY-EW effects 
and combined them with the SUSY-QCD ones. In this way we have obtained the simulta- 
neous MSSM predictions for both cr(pp(gg) — tc-|-tc) and B{h — S7) in full compatibility 
with the experimental data on the low-energy process b 37. The inclusion of the latter 
is crucial since it carries essential information on the FCNC parameters in the b-squark 
sector. This information can then be transferred to the LL part of the t-squark sector via 
CKM rotation and it thus enters in full interplay with the study of the production process 
PP(gg) ^ tc -I- tc. 

The aforesaid low and high energy interplay is particularly important in the context 
of the MSSM, where there are different and powerful sources of FCNC effects potentially 
inflecting the final prediction of B(h 37) [58,61]. For example, in the electroweak sector, 
we have the contributions from charged Higgs bosons and charginos, which have opposite 
signs and can have similar impact on the b ^ S7 amplitude. Although the charged Higgs 
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bosons turn out to contribute very little to pp(gg) tc + tc, their effect on b ^ S7 can be 
very important and, therefore, must be well balanced against the chargino contributions 
which, in contrast, affect the process pp(gg) — tc+tc in a significant manner. Only through 
a careful balance of the SM contribution against the various (strong and electroweak) 
supersymmetric effects (which may carry different signs) is possible to have a spectrum of 
relatively light SUSY particles and, at the same time, a safe value of that branching ratio 
within the experimental bounds. The ensuing supersymmetric spectrum obtained in this 
way, with some of the particles being relatively light, can then play a fundamental role to 
powerfully enhance the single top quark production process pp(gg) — tc + tc at the LHC. 

Such enhancement has to be measured with respect to the corresponding standard 
model (SM) vahic. To this end, we have computed the cross-section for pp(gg) — ^ tc + tc 
in the SM and found that it is extremely small (~ 10^^ pb). Since the supersymmetric 
effects can increase this result up to 7 orders of magnitude, it follows that the process under 
consideration could have a cross-section of order of 1 pb within the MSSM. This situation 
would be attained in the most favorable scenario, viz. in the presence of relatively light 
gluinos of a few hundred GeV and large flavor-mixing parameters of order S23 ~ 0.5 — 1.0. 
We note that a cross-section value of 1 pb amounts to 10^ events per 100 fb"^ of integrated 
luminosity. 

We should recall that the FCNC Higgs boson decays within the MSSM {H^,A^ 
tc,tc) can also be an additional source of single top quark events. However, as revealed by 
the analysis of [32], these FCNC decay modes are not competitive with the direct FCNC 
production process that we have studied here, the decay mechanism giving rates some 
two orders of magnitude smaller under the same set of assumptions. Furthermore, non- 
supersymmetric extensions of the Higgs sector of the SM, such as the general Two-Higgs 
Doublet Model (2HDM), do not provide any direct production channel at one-loop. In this 
context, the only possible source of tc and tc events comes from Higgs-boson decays, which 
may trigger FCNC decay modes owing to the enhanced charged Higgs boson interactions. 
These have been studied in [29, 30] and the result is that this source of single top quark 
final states is, again, not competitive with the direct mechanism pp(gg) ^ tc + tc within 
the MSSM. The latter, therefore, could be the most efficient source of single top quark 
FCNC events in renormalizable perturbartive extensions of the SM^. 

Apart from extending the results of [43] to include the full plethora of loop corrections 
originating from the SUSY-EW sector, in this work we have improved the calculation of 
the SUSY-QCD part in several respects as compared to the previous reference, mainly 
through the use of the rcnormalization group running values of the SM parameters, and 
also through an up-to-date set of PDF functions. Equipped with these updatings, we 
have compared the SUSY-QCD and SUSY-EW parts, separately and also in combination, 
and have identified and evaluated the region of the parameter space where the SUSY-EW 
effects can be sizeable in full compliance with the restrictions imposed by the low-energy 
observable B{h — > S7). That region, which is characterized by relatively light charginos 
and neutralinos, could be important, especially if the SUSY-QCD effects turn out to be 

^Othcr extensions of the SM, such as topcolor models, could give larger amount of tc pairs at the LHC, 
see e.g. Refs. [69,70]. 
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suppressed - e.g. if the gluinos are very heavy, say, of order of a few TeV. In this kind of 
scenarios, obviously controlled by the maximal SUSY-EW contributions, we find values of 
cTtc substantially smaller than those obtained under optimal SUSY-QCD conditions, but 
the event rates can still be sizeable (see below). 

The following observation is in order. The FCNC effects in the MSSM may involve, 
in general, all kinds of chirality mixings (LL, LR and RR) in the squark mass matrices, 
and all of them could be important. This was explicitly shown at the level of FCNC top 
quark decays in [19], and confirmed also for top quark production in [45,46]. However, 
in this work, following and extending the results of [43], we have limited ourselves to 
consider flavor mixing in the LL sector only. This is perhaps the most conservative point 
of view that one can adopt and, moreover, is well-motivated by renormalization groups 
arguments. That restriction notwithstanding, it is remarkable that the cross-section for 
pp(gg) tc + tc, within the SUSY-QCD favored scenario, can reach the ~pb level at the 
LHC in full compatibility with the experimental limits on the radiative B-meson decay. 
This is one of the most important results that we have obtained. 

Another important result is that, if the gluinos would be very heavy (say, with masses 
of order of a few TeV) and/or the misalignment of the quark and squark mass matrices 
would be negligible (i.e. 623 ~ 0), the total MSSM contribution could still be sizeable 
thanks to the potential SUSY-EW effects. In these conditions, the MSSM single top 
quark production cross-section could reach ~ 0.01 pb. This result, although two orders of 
magnitude smaller than the optimal SUSY-QCD one, entails 10'^ FCNC top quark events 
per 100 fb"^of integrated luminosity at the LHC and, therefore, it is still sizeable. 

Finally, we have also estimated the possible impact on the FCNC cross-section under 
study if we would introduce (on top of the b ^ 37 constraints) the restrictions induced by 
the experimental data on — B^ oscillations. We have found that it could be significant, 
namely it could lead to a reduction by roughly two orders of magnitude of the maximum 
result, hence a cross-section of order lO^^pb at most. However, we emphasize that the 
calculation of the SUSY contributions to the B^ — B^ mixing existing in the literature [49] 
has been obtained within the mass insertion approximation. They are, thus, only valid for 
sufficiently small values of 62^, namely for values quite away from our preferred value of 
this parameter, see Eq. ([4.71) . As a result, our optimized FCNC cross-section (~ Ipb) need 
not necessarily undergo the level of reduction mentioned above. The situation could only be 
reliably assessed if we would perform a computation of the SUSY contribution to B^ — B^ 
beyond the mass insertion approximation. But this is not the main aim of our work and, 
moreover, there are additional, intrinsic, sources of uncertainty in this calculation which are 
associated to our insufficient knowledge of the strongly interacting matrix element involved 
in the B^ — B^ amplitude. These irreducible uncertainties would remain unsettled. For 
this reason we have clearly separated the analysis of the two kind of constraints imposed 
on our calculation respectively by the low-energy processes b — > S7 and B° - B^, in the 
sense that they involve different sort of assumptions, those from B^ — B^ being susceptible 
of a higher level of uncertainty. 

Taking into account the important signature carried by the top quark and the virtual 
absence of FCNC contributions both from the SM and the general 2HDM, the detection 
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and appropriate identification of a buncli of events of this sort could be interpreted as a 
hint of SUSY physics. Let us note that the detection of tc-pairs in the LHC at the level 
< Ipb will be difficult [39,71], and that further experimental studies (including compar- 
ison of momentum distributions, etc.) will be necessary to assess the feasibility of such 
measurement. Such studies are beyond the scope of the present work. 

To summarize, SUSY interactions could be an efficient source of tc and tc events emer- 
ging from direct production FCNC processes at the LHC. These events, if effectively tagged 
and confidently discriminated among the other sources of single top quark production, 
could hardly be attributed to any other alternative FCNC mechanism within renormaliza- 
ble quantum field theory. Further studies will be necessary to fully assess this interesting 
possibility from the practical point of view, but in the meanwhile we have shown here 
that, if we use fairly conservative assumptions on the theoretical sources of flavor mixing 
in the MSSM and consider also the important phenomenological restrictions imposed by 
low-energy FCNC processes, there is still a potentially large supersymmetric enhancement 
of the FCNC single top quark signal at the LHC. 
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